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ABSTRACT 

A comparative study of raindrop size distribution measurements has been conducted at NASA’s Goddard 
Space Flight Center where the focus was to evaluate the performance of the upgraded laser-optical OTT 
Particle Size Velocity (ParsiveP; P2) disdrometer. The experimental setup included a collocated pair of 
tipping-bucket rain gauges, OTT Parsivel (PI) and P2 disdrometers, and Joss-Waldvogel (JW) disdrometers. 
Excellent agreement between the two collocated rain gauges enabled their use as a relative reference for 
event rain totals. A comparison of event total showed that the P2 had a 6% absolute bias with respect to the 
reference gauges, considerably lower than the PI and JW disdrometers. Good agreement was also evident 
between the JW and P2 in hourly raindrop spectra for drop diameters between 0.5 and 4 mm. The P2 drop 
concentrations mostly increased toward small sizes, and the peak concentrations were mostly observed in the 
first three measurable size bins. The PI, on the other hand, underestimated small drops and overestimated the 
large drops, particularly in heavy rain rates. From the analysis performed, it appears that the P2 is an im- 
provement over the PI model for both drop size and rainfall measurements. P2 mean fall velocities follow 
accepted terminal fall speed relationships at drop sizes less than 1 mm. As a caveat, the P2 had approximately 
1 m s^^ slower mean fall speed with respect to the terminal fall speed near 1 mm, and the difference between 
the mean measured and terminal fall speeds reduced with increasing drop size. This caveat was recognized as 
a software bug by the manufacturer and is currently being investigated. 


1. Introduction 

The Particle Size Velocity (Parsivel) is a laser-optical 
disdrometer that measures the size and fall velocity of 
hydrometeors. Ldffler-Mang and Joss (2000, hereafter 
LJOO), who introduced the specifications and measure- 
ment principles of the Parsivel, described three major 
advantages with respect to commercially available dis- 
drometers at that time. First, they pointed out the easy- 
to-operate, robust, and low-cost features of the Parsivel. 
Second, they mentioned that the Parsivel could estimate 
the size of drizzle drops down to a 0.1-mm diameter with 
modifications made to the optical system. Third, they 
stated that the Parsivel estimates the size and fall ve- 
locity of snowflakes and is useful for discriminating the 
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hydrometeor type, which makes the Parsivel a present 
weather sensor. 

LJOO evaluated the performance of the Parsivel 
through comparison with a collocated impact-type 
Joss-Waldvogel (JW) disdrometer (Joss and Waldvogel 
1967) and a recording Hellmann accumulation gauge 
(Sevruk 1996) in rain. They found reasonable agreement 
between the JW and Parsivel raindrop size distribution 
(DSD) in the 0.7-2-mm drop-diameter range, while the 
Parsivel detected higher concentrations at sizes less 
0.7 mm, and a limited sample was the cause for the no- 
ticeable differences for drops larger than 2 mm in diameter. 
Reasonable agreement was also noted between the 
JW and Parsivel rain-rate time series as well as between 
Parsivel and the Hellmann gauge for the daily rain totals. 

The Parsivel disdrometer became commercially avail- 
able through PM Tech Inc., and Yuter et al. (2006) op- 
erated a loan unit in a rain/wet snow event and in a dry 
snow event in the Cascade and Rocky Mountains, 
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respectively. They found that the hydrometeor fall ve- 
locities were noticeably different at different tempera- 
ture ranges as the phase of hydrometeors change. The 
size distributions were also noticeably different during 
different storm phases as classified by environmental 
temperature. 

Two PM Tech Parsivel units were operated at Na- 
tional Aeronautics and Space Administration (NASA)’s 
Wallops Flight Facility (WFF) for several years, where 
both units recorded 12% and 19% less rainfall than 
collocated rain gauges in 30 rain events. These two units 
also participated in the Canadian CloudSatl Cloud- 
Aerosol Lidar and Infrared Pathfinder Satellite Observa- 
tions (CALIPSO) Validation Project (C3VP). Battaglia 
et al. (2010) used data from C3VP to evaluate the 
performance of the PM Tech Parsivel in snow relative 
to a collocated two-dimensional video disdrometer 
(2DVD). They found significant differences in snow- 
flake size distribution and fall velocities between the 
two instruments and focused on the shortcoming of 
Parsivel measurements. 

Krajewski et al. (2006) compared a PM Tech Parsivel 
with a dual-beam spectropluviometer (DBS) (Hauser 
et al. 1984) and a 2DVD (Kruger and Krajewski 2002; 
Schonhuber et al. 2007) in rain. Agreement between the 
DSD of the three disdrometers was observed for drop 
diameters ranging from 0.5 to 4 mm, while the Parsivel 
had lower drop concentrations at sizes less than 0.5 mm 
in five major events. The mean Parsivel fall velocity was 
higher than terminal fall speeds for raindrops (Beard 
1976) at sizes less than 1 mm in diameter, and the reverse 
was true at sizes above 1.4 mm. Krajewski et al. (2006) 
reported that the Parsivel estimated higher rain accu- 
mulations than both of the other collocated disdrometers 
as well as nearby tipping-bucket rain gauges in each of 
the five major events, but the differences were 16% or 
less with respect to the rain gauges. 

In 2005, OTT Hydromet purchased all the rights of 
the Parsivel and redesigned the instrument to specifi- 
cally target applications in the transportation sector. 
Indeed, transportation agencies around the world pur- 
chased the OTT Parsivels (PI) to be used as a present 
weather sensor. There are several differences between 
the OTT and PM Tech models. The output voltages in 
the PM Tech and OTT Parsivels were sampled at 10 and 
50 kHz, respectively (Battaglia et al. 2010). One sample 
is therefore taken every 0.1 and 0.02 ms in the PM Tech 
and OTT Parsivels, respectively. The PI used a much 
less expensive laser device and reduced homogeneity 
across the laser beam resulted in an overestimation of 
large drops relative to a collocated 2DVD in northern 
Alabama (Tokay et al. 2013, hereafter TPGW13). The 
event rain totals from PI were nearly 14% higher than 


estimates from collocated tipping-bucket gauges. The 
PI was calibrated using a rotating precision disk, and the 
mean fall velocities at sizes above 1.4 mm in diameter 
were lower but much closer to the terminal fall speed 
than the PM Tech Parsivel. 

The PI participated in the first field experiment for 
intercomparison of rain intensity measurements; the 
field experiment was organized by the World Meteoro- 
logical Organization in central Italy, where 30 selected 
rain gauges and optical disdrometers operated for 
1.5 years (Vuerich et al. 2009). The PI rated a score of 3 out 
5 for 1-min rain rate with respect to selected reference 
gauges. Lanza and Vuerich (2009) pointed out that op- 
tical disdrometers overestimated rain intensity and the 
overestimation increased with increasing rain rate. This 
finding is consistent with the study by Thurai et al. (2011) 
and TPGW13. 

In 2011, OTT introduced a new model of the Parsivel, 
namely, Parsivel^ (P2). A beta version of P2 was eval- 
uated through a comparative study at NASA’s Goddard 
Space Flight Center. This study presents the perfor- 
mance of P2 with respect to the PI model through an 
experimental study in rain. The experimental study also 
includes a JW disdrometer and two tipping-bucket rain 
gauges. It should be noted that none of the disdrometers 
used in this study is considered as an absolute reference. 
The technical details of the Parsivel disdrometer are 
given in section 2, which is followed by the rainfall sta- 
tistics and the measurement site and data in sections 3 
and 4, respectively. The comparison of event rainfall 
totals between the disdrometers and rain gauges is also 
given in section 4. The comparison of hourly DSD be- 
tween the old and new models of the OTT Parsivel and 
collocated JW is presented in section 5, while section 6 
provides the comparison of fall velocity measurements 
between the old and new OTT Parsivel. The probability 
of and cumulative distributions of fitted gamma pa- 
rameters and integral rain parameters are presented in 
section 7 followed by conclusions in the last section. 

2. Parsivel disdrometer 

All the Parsivel models measure the size and fall ve- 
locity of hydrometeors that fall through a laser sheet 
independently. The laser sheet is approximately 180 mm 
long, 30 mm wide, and 1 mm high. The size of the hy- 
drometeor is estimated from the maximum attenuation of 
the signal. The duration of the hydrometeor within the 
laser beam provides the fall speed. The raindrops falling 
through the laser sheet are assumed to be spherical for 
sizes less than 1 mm in diameter, while the drop axis ratio 
varies linearly from 1 to 0.7 for drops between 1 and 
5 mm. For drops larger than 5 mm, the axis ratio is set to 


1278 


JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY 


Volume 31 


0.7. The Parsivel processing software assumes snowflakes 
as spheres and therefore the measured “Parsivel size” is 
a one-dimensional length, which is not necessarily rep- 
resentative of the equivalent diameter. 

The Parsivel’s nominal sampling area is 54 cm^, but 
the effective sampling area is calculated considering 
the partially observed hydrometeors. While the sampling 
area was intentionally designed to reduce the probability 
of multiple hydrometeors being detected at the same 
time, this may not be avoidable, particularly in the pres- 
ence of abundant small hydrometeors. Also, secondary 
particles may result from splashes, and windy conditions 
can result in the particles falling through the laser sheet at 
an angle. For analysis of raindrops herein, we ehminated 
the suspect false particles whose fall speeds fell outside 
the ±50% of the generally accepted range of drop ter- 
minal fall speed; however, this subjective criterion elim- 
inates mainly the small drops. The terminal fall speed 
table that is used in this study is based on wind tunnel 
experiments following Beard (1976). 

The Parsivel’s raw output provides the number of 
drops in a 32 X 32 size versus fall velocity matrix. The 
size range is from 0 to 25 mm, and the class width in- 
creases with the size from 0.125 to 3 mm. The first two 
size classes are left empty due to the low signal-to-noise 
ratio, and the minimum detectable size is approxi- 
mately 0.25 mm. The fall velocity range is from 0 to 
20 m s~^ and the class width increases with fall velocity. 
Thus, the Parsivel cannot provide precise size and fall 
velocity of the hydrometeors, particularly at larger sizes 
and faster fall speeds. The calculated DSD parameters 
and maximum size, therefore, are subject to quantiza- 
tion errors. 

Homogeneity of the laser sheet is crucial for accuracy 
of the Parsivel measurements. LJOO conducted a labo- 
ratory study to test the accuracy of size and fall velocity 
of PM Tech Parsivel measurements. They showed that 
the overall error in estimating the drop diameter does 
not exceed ±0.1 mm plus ±5%, and the errors in fall 
velocity are within 25% for 0.3-mm-diameter drops and 
10% for 5-mm-diameter drops. High-precision refer- 
ence particles at 0.5, 1, 2, and 4 mm are employed to 
calibrate PI with a rotating disk. The manufacturer re- 
ports that the accuracy of rain rate is ±5%, while the 
radar reflectivity accuracy is ±20%. A key issue with PI 
was the inhomogeneous laser beam. It was pointed out 
that large particles could be overestimated by 20%, 
depending on the precipitation intensity (K. Nemeth, 
OTT, 2008, personal communication). Indeed, Thurai 
et al. (2011) showed higher median mass diameters 
(Dmass) estimated by the PI than the 2DVD at D^ass > 
2 mm, and when the rain intensity was greater than 
20mmh“^. TPGW13 also showed that higher raindrop 


concentrations were detected by the PI compared to the 
JW and the 2DVD. The PI raindrop fall velocities, on 
the other hand, agreed better with the terminal fall 
speed of raindrops, particularly for midsize drops (1- 
3 mm) than those from PM Tech (TPGW13). The sam- 
pling output rate is user selectable but can be as high as 
10 s as opposed to 1-min output sampling for the PM 
Tech Parsivel. 

P2 is a new design where the electronics are placed in 
the sleeve rather than in the sensor’s two heads as with 
the PI. The sensor head is also heavier than the previous 
model’s. The manufacturer suggests that the key im- 
provement is the use of a more expensive laser device 
and better homogeneity of the laser sheet. The mea- 
surement accuracy was reported to ±1 size class in P2 as 
opposed to ±3 class in PI up to 2 mm, while for the 
particles above 2 mm, the accuracy is ±0.5 size class in 
P2 as opposed to ±2 in PI. It should be noted that the 
class width doubles at 2.8 mm. 

3. Rainfall statistics 

Percent bias and percent absolute bias are used to 
compare the event rain totals between the disdrometer 
and rain gauges, and raindrop size distribution and in- 
tegral rainfall parameters between the disdrometers. 
The percent bias and absolute percent bias between the 
measurements of the two instruments (x, y) for n sam- 
ples are calculated as 


, . bias 

percent bias = -^^ u, 

{x,y) 

(1) 

absolute bias 
percent absolute bias = , ^ , 

{x,y} 

(2) 

where 


1 " 

bias = - X (a- 
n ,=i 

(3) 

1 " 

absolute bias = - ^ x.— y. , and 
n ‘ 

(4) 

(^a) L \ ■ 

n 2 

(5a) 

If one of the variables is from a reference instrument 
(e.g., x), then the denominator in Eqs. (1) and (2) is 

expressed as 


1 " 

w= - 

(5b) 
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Fig. 1. In situ precipitation measuring test site at the roof of 
building 33 at NASA’s Goddard Space Flight Center, Greenbelt, 
MD. OTT Parsivel and ParsiveP are on the post at the left and 
right sides, respectively. The impact-type JW disdrometer is in the 
middle sitting on a cinder block with green water-absorbing ma- 
terial. Two tipping-bucket rain gauges are in the middle behind the 
JW disdrometer sitting on a wooden box. A third tipping-bucket 
gauge was also in the field (left) but did not operate for the ex- 
periment period. 

In this study, the gauge that was operated for all rain 
events is considered as a reference. The tipping-bucket 
rain gauges are manufactured by Met One Inc. (model 
370) and were previously tested through comparative 
field studies with Texas Electronics (model TR-525I) 
rain gauges in northern Alabama (TPGW13), and with 
Hydrological Services (model TB3) syphon rain gauges 
at NASA WFF, Wallops Island, Virginia. 

4. Measurement site and data 

During tests of the P2, the P2, PI, JW, and two Met 
One Inc. tipping-bucket rain gauges were collocated on 
the roof of building 33 at NASA’s Goddard Space Flight 
Center (39.00°N, 76.87°W) (Fig. 1). The time of each tip 
(0.254-mm accumulation) was recorded for the gauges, 
while drop counts detected by all three disdrometers 
were first recorded at 10-s intervals and then integrated 
to 1-min values. The impact-type JW disdrometer de- 
termines the raindrop size by assuming that the rain- 
drops are falling at their respective terminal fall speeds 
and provides the drop counts in both 127 and 20 unequal 
size classes. We used the former to calculate the integral 
rain parameters and the latter to compare the composite 


DSD with the Parsivel disdrometers. The details of JW 
disdrometers can be found in Tokay et al. (2005) and the 
references therein. 

During the experiment period (28 February-2 August 
2011), there were 36 rain events where the rain totals 
exceeded 1-mm rainfall. Here, we define a rain event as 
a rain period separated by 2-h or longer rain-free pe- 
riods in the rain-rate time series of the JW disdrometer, 
and a rain/no-rain threshold was set at a minimum of 10 
drops and a rain rate of O.lmmh”^ for all three dis- 
drometers. For total rainfall, the reference rain gauge 
recorded 381 mm of rainfall. The second gauge failed 
to operate in four events; nonetheless, it recorded just 
2% less rainfall than the reference gauge when the 
total rainfall was recalculated from the remaining rain 
events. In comparisons to the gauges, the P2 recorded 
349 mm of rainfall in 145 rainy hours, while PI and JW 
had 404 mm in 134 rainy hours and 325 mm in 166 h, 
respectively. 

Considering event rain totals, excellent agreement 
was found between the two rain gauges, where percent 
bias and percent absolute bias were —2% and 4%, re- 
spectively, in 32 events (Fig. 2a). The JW disdrometer 
recorded higher event rain totals than the rain gauge 
when the event total was less than 4 mm and the con- 
verse was true at higher event totals (Fig. 2b). The ab- 
solute bias was 15%, higher than a similar study in 
northern Alabama (TPGW13), but it was within the 
manufacturer-specified limits (Tokay et al. 2005). The 
PI overestimated half of the rain events with respect to 
the gauge and the percent absolute bias was 18%, higher 
than the TPGW13 study (Fig. 2c). The P2 had the best 
agreement with the rain gauge in event rain totals, where 
the absolute bias was 6% (Fig. 2d). The differences in 
rainfall characteristics and environmental conditions 
resulted in differences in rainfall statistics between this 
and the TPGW13 study. It seems possible that wind and 
turbulent conditions on the roof could have contributed 
to higher differences in event totals between the in- 
struments. Unfortunately, there was no recording wind 
sensor available. However, wind measurements from 
a nearby station are provided for the comparison of the 
DSD measurements in the next section. 

The Parsivel and JW disdrometers have nearly the 
same sampling area and we used the same rain/no-rain 
criteria, so the comparison of rainy minutes between the 
disdrometers provides information on the sensitivity of 
each disdrometer to rain. The JW disdrometer recorded 
more rainy minutes than both Parsivels in every rain 
event, and the percent absolute bias was 20% and 12% 
between JW and PI and between JW and P2, respec- 
tively (Figs. 3a,b). The P2 recorded more rainy minutes 
than PI, where the absolute bias was 8% (Fig. 3c). This 
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Fig. 2. Comparison ol event rain totals (a) between the two MetOne Inc. tipping-bucket rain gauge, and (b) between the reference 
gauge and JW, (c) between the reference gauge and the OTT Parsivel, and (d) between the reference gauge and OTT ParsiveP dis- 
drometers. The percent bias and absolute bias are also given. 


demonstrates that P2 is more sensitive than PI but not as Considering the Parsivel disdrometer, the raw output is 
sensitive as the JW disdrometer for detecting rainfall. the number of drops at the /th size and /th velocity bin 

(C/y) and DSD is calculated as 


5. Raindrop size distribution measurements 

Midsize (1-3 mm) drops are the main contributor to 
rainfall, while large drops (greater than 3 mm) contrib- 
ute significantly in heavy rain. The agreement between 
the disdrometer DSDs at these size ranges is directly 
linked to the differences in event rain totals discussed 
previously. The DSD is expressed as the number of 
drops per volume of air for a given drop size interval. 


N(D) 


1 

Time 


c.. 


Ji^?iu.Area(D,)AD,’ 


( 6 ) 


where AD, is the width of the ith size bin; and n and m 
are the number of size and velocity bins, respectively, 
and both are equal to 32. Area (D,), the effective dis- 
drometer sampling area, is calculated considering par- 
tially detected drops across Parsivels’ laser sheet and is 



Rainy minutes 


Rainy minutes 


Rainy minutes 


Fig. 3. Comparison of event rainy minutes (a) between the JW and OTT Parsivel disdrometers, (b) between the JW and OTT ParsiveP 
disdrometers, and (c) between the OTT Parsivel and OTT ParsiveP disdrometers. The bias and absolute bias are also given. 
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Fig. 4. Comparison of hourly raindrop spectra from 20 different rain events. The hourly spectrum is given for the JW disdrometer 
(solid), the OTT Parsivel disdrometer (dashed), and the OTT ParsiveP disdrometer (dashed-dotted). The hourly mean and maximum 
wind speeds are also given (kmh^'). 


equal to 180 mm X (30 mm — Z),/2). The observational 
period, Time, is set to 60 s in this study. Although vj is the 
measured raindrop fall speed at the jth velocity bin in 
Eq. (6), the terminal fall speed u(D,) is used in this study 
to eliminate any instrument-measured fall speed arti- 
facts on the DSD, and therefore it provides a more direct 
comparison with JW-derived DSD. To avoid sampling 
fluctuations, we make sure that all composite DSDs 
have more than 10 rainy minutes in an hour and that the 


differences in the disdrometer’s rainy minutes are less 
than 15%. Figure 4 shows 30 one-hourly DSD compos- 
ites for the JW, PI, and P2 for 20 rain events. 

One shortcoming of the PI is the underestimation of 
small drops with a sharp drop-off occurring toward 
smaller sizes after peak concentrations at 0.84 or 
0.96 mm in diameter (TPGW13). In this study, PI had 
peak concentrations at the same drop diameters except 
for a few rainy hours where the peak concentration 
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occurred at a diameter of 0.58 mm. The sharp drop-off 
toward smaller sizes was mostly observed in PI hourly 
composites. The P2 exhibited peak concentrations mostly 
between 0.32 and 0.58 mm, with no sharp decrease to- 
ward the smaller sizes. The JW disdrometer, on the other 
hand, had peak concentrations mostly from 0.44 to 
0.65 mm but showed drastic drop-offs toward small sizes 
in a few cases (Figs. 4m; a, [B, y). 

For midsize drops, the P2 and JW exhibited good 
agreement, while the PI indicated higher concentrations 
of drops larger than 2.44 (Fig. 4; a) and 2.82 mm (Fig. 4; 
0). The overestimation of mid- to large-size drops has 
been previously reported as another shortcoming of PI 
(TPGW13). The P2 also exhibited higher concentrations 
of larger drops than the JW in a number of hourly 
composites (Figs. 4k,l; y), but the differences in con- 
centration were relatively less between the JW and P2 
compared to the JW and PI. Although we did not select 
any disdrometer as a reference, the good agreement 
between the JW and P2 hourly DSDs in a broad size 
range (0.5-4 mm) suggests that P2 is likely a better in- 
strument in measuring the size spectra than PI. 

For the size spectra at drop diameters larger than 
4 mm, the sampling is limited and, of course, the JW 
cannot distinguish drops larger than 5.3 mm. This makes 
it difficult to evaluate the full performance of the P2. 
However, the slope of the P2 DSD often followed JW, 
while PI measured distinctly higher concentrations at 
these sizes (Figs. 4w,x; a, fi). It should be noted that for 
several hourly DSDs, all three disdrometers showed 
good agreement except at the lower end of the spectrum 
in a number of hourly DSDs (Figs. 4b-d,q,s-v). For 
hourly DSDs, where the differences were significant, 
there was no wind dependency as previously suggested 
in TPGW13. 

6. Fall velocity measurements 

Raindrop trajectories may deviate from the vertical in 
windy conditions, but it is expected that the raindrop 
mean fall speed should follow the terminal fall velocity 
quite closely. At least the mean fall speed should not 
systematically deviate from the terminal fall velocity for 
all rain events. Figure 5 shows the mean and standard 
deviation of PI and P2 fall velocities for 12 rain events, 
which recorded at least 10 mm of rainfall. The terminal 
fall velocity is also depicted as a reference. 

There were significant differences in fall velocity 
measurements between the two disdrometers. PI over- 
estimated the fall velocities of small drops, where the 
difference between the mean measured and terminal 
fall velocity peaked at the smallest measurable size 
(0.32 mm) and the difference from terminal fall velocity 


was as high as 42% (Fig. 5a). The overestimation di- 
minished with increasing size and mean fall velocities 
were less than terminal fall velocities, where the cross- 
over occurred between 0.95 (Fig. 5k) and 1.67 mm (Fig. 
5g). The difference between the mean measured and 
terminal fall velocity was as high as 18%, which occurred 
at 1.67 mm (Fig. 5k). This was also a fairly windy event, 
with a maximum wind speed of approximately 32 km h“ 
In a number of events, PI overestimated the fall velocity 
of large drops at sizes 3.35 mm or larger (Fig. 5w). The 
difference between the mean measured and terminal fall 
velocity was as high as 32%, which was occurred at 
6.69 mm (Fig. 5w). 

The P2 mean fall velocities followed the terminal fall 
speed at smaller drop diameters up to 1.09 mm, where 
differences between measured and terminal fall speeds 
remained less than 15%. At diameters of 1.09mm, the 
P2 underestimated the fall speeds by 12% (Fig. 5j) and 
27% (Fig. 5n). The manufacturer recognized this short- 
coming as a software bug and is currently investigating 
the issue (K. Nemeth, OTT, 2013, personal communica- 
tion). Underestimation of fall velocities was also evident 
for larger drops, but the difference between the measured 
and terminal fall speeds was not as pronounced. In a 
number of events, P2 overestimated the fall velocity of 
the drops at sizes larger than 3.35 mm. 

Fall velocity measurements play an important role in 
calculating the DSD, as shown in Eq. (6), as well as in 
rain parameters, except in rain rate. Figure 6 shows the 
6-hourly P2 DSD from six different events where the 
DSD was calculated from both measured and terminal 
fall speeds. Visually, the differences between the mea- 
sured and terminal fall speed DSDs do not appear sig- 
nificant; however, a closer look at the 1.09-mm bin 
revealed that concentrations were 23% lower when the 
terminal fall speeds were used. This is consistent with 
approximately lms“^ differences between the mean 
measured and terminal fall speeds. 

Table 1 presents statistics that quantify the role of fall 
speed in total concentration, Nt\ liquid water content, 
W\ reflectivity; and mean The statistics are based 

on 8685 P2 observations, each of 1 min, from the 36 rain 
events mentioned previously. The underestimation of 
fall velocities in the midsize drop-diameter bins domi- 
nated the DSD and rain parameters, resulting in a nega- 
tive percent bias, since these parameters were calculated 
using terminal fall speed first. As in Eq. (6), the fall speed 
is in the denominator in calculating the DSD and integral 
rain parameters and, therefore, the underestimation of 
measured fall speeds results in higher DSD and rain 
parameters. The percent absolute biases for Nt and W 
were 8% and 13%, respectively, while Z had a 0.6-dB 
absolute bias. The mean mass diameter the ratio 
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Fig. 5. Mean (squares) and standard deviation (vertical bars) of measured fall velocity by the 
OTT Parsivel and the ParsiveP disdrometers as a function of diameter for 12 rain events. The 
terminal fall speed following Beard (1976) is shown as a solid line. The event mean and 
maximum wind speeds are also given (kmh^^). 
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Fig. 6. Hourly OTT ParsiveP raindrop size distributions from six different events. Hourly size distributions are calculated using measured 

(dashed) and terminal (solid) fall speeds. 


of the fourth to the third moment of the DSD, exhibited 
a relatively insignificant bias of 0.02 mm. 


N{D) = NJJm) 


D 


D„ 


exp 


- (4 + m) 


D 


D„ 


, ( 8 ) 


7. Raindrop size distribution and integral rainfall 
parameters 

Probability and cumulative distributions of DSD and 
integral rain parameters were calculated from 8459 ob- 
servations, each of Imin, for which all three dis- 
drometers reported rainfall. To be consistent with the 
JW disdrometer, the terminal fall speed was used in cal- 
culating DSD and rainfall parameters from the Parsivel 
disdrometers. The DSD is formulated by a three-parameter 


IS given as 


N{D)=N*f^im) 


D 




exp 


— (4 + m) 


D 




where N* and Nw are the normalized intercept pa- 
rameters with respect to Nt and W, respectively; 
and m is the shape parameter. The terms fi and fZ 


Table 1. Rainfall statistics of raindrop size distribution 
(Dmass) and integral rain {Nt, W, Z) parameters that are calcu- 
lated from OTT ParsiveP measured and terminal fall speeds. 
Bias and percent bias are calculated using terminal fall speed 
first. 




Bias 

ab_bias 

per_bias (%) 

per_ab_bias (%) 


Nt (m ^) 

-26 

35 

-6 

8 

and 

W (gm^^) 

-0.017 

0.625 

-13 

13 


Z(dB) 

-0.6 

0.6 

-1 

2 

(7) 

Hmass (mm) 

-0.01 

0.02 

-1 

2 
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are the function of the shape parameter and are ex- 
pressed as 


, (4 + m)'"+l 

r(™ + i) ““ 

(9) 

6 (4 + mr^ 

^2'' ^ 256 r(m + 4) ’ 

(10) 


respectively. Both normalized intercept parameters and 
Dmass were directly calculated from disdrometer mea- 
surements, while the shape parameters were obtained by 
minimizing the error between the observed and gamma- 
fitted DSD-based rain rates. The complete gamma 
function was assumed in calculating the rain rate from 
the gamma-fitted DSD. The range of the shape param- 
eter was considered to be between —2 and 20, even 
though higher values were occasionally reported (Tokay 
and Short 1996). Approximately 12% of the data were 
disregarded because they were outside this range. Out- 
of-range values for the shape parameters are partly due 
to the use of complete gamma functions and partly be- 
cause the gamma fit was not well suited to the observed 
DSD spectra. 

Probability and cumulative distributions of integral 
parameters {R, W, Z) and D^^ss showed good agreement 
between the three disdrometers (Figs. 7a-h), but there 
were noticeable differences in mean values of these 
parameters between the instruments (Table 2). This is 
the result of the differences in the disdrometer DSDs, 
where a particular size regime contributes the most to 
a particular rain parameter. The contributions of small, 
midsize, and large drops differed from each other for R, 
W, and Z at different rain intensities. The P2 had higher 
mean rain rates than JW and lower mean rain rates than 
PI. For the liquid water content, both Parsivels had 
similar means, which were higher than the JW mean 
value. P2 had 1 and 3.4 dB lower mean reflectivity than 
the JW and PI mean reflectivity, respectively. The un- 
derestimation of small to midsize drops resulted in lower 
mean W and R in JW disdrometers (Figs. 4m; a, y) and 
the overestimation of midsize to large drops resulted in 
higher R and Z in PI (Figs. 4k; a, (3). 

Probability and cumulative distributions of intercept 
and shape parameters of the gamma DSD showed good 
agreement between JW and P2, but the distributions 
were shifted toward lower intercept and higher shape 
parameters in PI (Figs. 7i-p). P2 had the highest mean 
intercept parameter, while the PI had the lowest. This 
was due to the underestimation of small drops in the JW 
and PI, where the latter exhibited a sharp drop off at 
small sizes (Figs. 4j,m,o,w). The mean shape parameters 
of the P2 were 5.6 and 5.3 for fitted gamma DSDs based 


on Eqs. (6) and (7), respectively, in a good agreement 
with the mean shape parameter of the JW. 

The distribution of the triplet of gamma DSD pa- 
rameters plays an important role in the formulation 
of the spaceborne radar rainfall retrieval algorithms. 
NASA’s Tropical Rainfall Measuring Mission (TRMM) 
precipitation radar rainfall algorithm employs a gamma 
DSD with fixed shape parameter (m = 3) (Kozu et al. 
2009), also indicated in the probability distributions in 
Figs. 7k, 1. NASA’s Global Precipitation Measurement 
(GPM) mission precipitation radar rainfall algorithm, 
on the other hand, will also use the gamma-based DSD, 
and an empirical relationship between the parameters of 
gamma DSD, such that it can retrieve all parameters of 
DSD from its dual-frequency measurements, is highly 
desirable (R. Meneghini, NASA Goddard Space Flight 
Center, 2013, personal communication). 

8. Conclusions 

Affordability, easy installation and maintenance, and 
robustness of hardware in the field have resulted in 
worldwide use of the Parsivel disdrometer in many ap- 
plications of Earth sciences and beyond. However, there 
have been relatively few independent studies that have 
evaluated the performance of the OTT Parsivel through 
field studies. As noted in TPGW13, PI underestimated 
the concentrations of small drops less than 0.76 mm and 
overestimated the drop concentration for drop diameters 
larger than 2.4 mm, particularly in heavy rain. The fall 
velocity measurements made using PI seem to be better 
than the PM Tech version, as the mean fall speeds 
measured by the former followed well-established em- 
pirically derived terminal fall speed relationships more 
closely. 

Herein, we compared a collocated newer version of 
the OTT Parsivel (Parsivel^) to its predecessor, a JW 
disdrometer, and a collocated rain gauge pair. The OTT 
Parsivel^ showed better agreement with the reference 
rain gauge than the Parsivel and JW disdrometers. Good 
agreement was also evident in hourly DSDs between JW 
and Parsivel^, with drop diameters in the range of 0.5- 
4 mm. Parsivel^ drop concentrations generally increased 
toward smaller size bins, peaking in the first three 
measurable size bins of 0.34-0.58 mm. These features 
demonstrated that the Parsivel^ is indeed an improved 
version of the OTT Parsivel for the raindrop size and 
rainfall measurements. The Parsivel^-measured rain- 
drop fall velocities approximately followed the expected 
terminal fall speeds at the small end of the drop spec- 
trum, but there was approximately lms“^ difference 
between the mean measured and terminal fall speeds at 
1.09 mm. This discrepancy may exert a pronounced 
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Fig. 7. Probablity and cumulative distributions of (a),(e) rain rate; (b),(f) liquid water content; (c),(g) reflectivity; 
(d),(h) mean mass diameter; (i),(m) intercept of gamma-fitted distribution with respect to Nt and D^nass', (j),(n) intercept 
of gamma-fitted distribution with respect to W and Dmassi (k)>(o) shape parameter of gamma-fitted distribution with 
respect to AJ.; and (l),(p) shape parameter of gamma-fitted distribution with respect to N^, from JW (solid), OTT 
Parsivel (dashed), and OTT ParsiveP (dashed-dotted) disdrometers. The shape parameter of 3, which is used for the 
TRMM precipitation radar algorithm, is marked as a vertical dash line in (k) and (1). 
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Table 2. Rainfall statistics of raindrop size distribution (Dmass, and integral rain (R, W, Z) parameters. The terms nii and 

ni 2 refer to the shape parameters of the gamma distribution with respect to AJ. and Nw 



mean; std dev (JW) 

mean; std dev (Parsivel) 

mean; std dev (ParsiveP) 

R (mmh 

2.2; 5.2 

2.9; 9.3 

2.5; 6.1 

W(gm-3) 

0.120; 0.210 

0.139; 0.340 

0.140; 0.255 

Z(dB) 

35.7; 40.0 

38.1; 47.9 

34.7; 42.7 

^mass (mrn) 

1.16; 0.42 

1.28; 0.51 

1.16; 0.45 

VJ. (m^^mm^^) 

332; 535 

180; 211 

465; 791 

Nw (m^^mm^^) 

8317; 18486 

4637; 6444 

10088; 19337 

mi 

5.7; 3.8 

8.1; 4.9 

5.6; 4.2 

m2 

5.7; 4.3 

7.1; 4.9 

5.3; 4.5 


effect on computed DSD and integral rainfall parame- 
ters. Though the manufacturer will eventually correct 
the fall velocity issue in its software, the impact of the 
error can be mitigated for rain, as it is feasible to use 
terminal fall speed in calculating DSD and rainfall pa- 
rameters with modest error. However, the fall speed 
error may not be as easy to mitigate in snow, since the 
flake size versus terminal fall speed relation depends on 
habit and environmental conditions. The NASA GPM 
Cold-Season Precipitation Experiment in southern 
Ontario, Canada, and the World Meteorological Orga- 
nization lead the Solid Precipitation Intercomparison 
Experiment in Finland and Canada will provide an ex- 
cellent opportunity to evaluate the performance of the 
Parsivel^ in snow. 

As reported in LJOO study, the Parsivel could estimate 
the size of drizzle drops down to 0.1 mm in diameter with 
a modification in the optical system. This will provide 
a better estimate of total concentration and improve 
DSD parameterization. The raw output could be the 
specific information for each drop, including its time 
stamp, size, and fall velocity, rather than 32 size X 32 fall 
velocity matrix. The problematic issue of fall velocity 
measurements at 1 mm and larger sizes is expected to be 
resolved through a software update. This study is limited 
to rainfall, and similar comparative studies between 
the collocated disdrometers are needed to evaluate the 
performance of the OTT Parsivel^ in snow. It should be 
pointed out that communications between the manu- 
facturer and the users is key for diagnosing and cor- 
recting any issues of the instrument as well as for 
planning the next generation. 
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